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Summary  of  Accomplishments 

Simulations:  effects  of  interfacial  films  on  thermal  conductance 

•  Found  that  inserting  an  interfacial  film  at  ideal  argon-“heavy  argon”  interfaces  can  increase  Hbd  by  as  much  as 
23%.  Optimum  enhancement  occurs  at  low  temperatures  (10%  of  Tmeit),  very  thin  films  (<  10  unit  cells),  and, 
all  else  equal,  when  the  film  has  an  atomic  mass  equal  to  the  average  of  the  abutting  materials  [1,2]. 

•  Found  that  small  degrees  of  interfacial  mixing  can  actually  further  increase  /ibd-  In  the  best  case,  enhancement 
to  the  baseline  Hbd  increased  from  23%  (sharp  interface)  to  53%  (“medium”  mixing)  [2]. 

•  Developed  a  new  interatomic  potential,  the  angular-dependent  embedded  atom  method  (A-EAM),  to  describe 
multicomponent  systems  of  both  metallic  and  covalent  materials  [3].  Parameters  have  been  published  for  the 
binary  systems  Au/Si,  Au/Ge,  and  Al/Si.  Parameters  for  the  ternary  Au/Al/Si  system  are  to  be  published. 

Theory:  understanding  conductance  via  phonon  transport  theory 

•  Developed  and  tested  modifications  to  the  diffuse  mismatch  model  (DMM)  for  predicting  Hbd-  including  optical 
phonons  [4],  including  full  phonon  dispersions  [5],  breaking  detailed  balance  [6],  and  including  anharmonic 
(“inelastic”)  phonon  processes  [7].  Modifications  were  shown  to  improve  the  predictive  ability  of  the  DMM. 

•  Using  Green’s  functions,  determined  analytical  expressions  that  capture  the  effects  of  interfacial  bonding  and 
impurities  on  phonon  transmissivity.  Used  those  results  to  identify  cases  for  optimum  transmission  [8]. 

•  Applied  the  wavelet  transform  to  atomistic  simulations  to  reveal  phonon  transport  processes  [9].  At  interfaces, 
this  approach  reveals  the  effects  of  mixing  and  roughness  on  phonon  scattering  dynamics. 

•  Using  the  interfering  particle  model,  analyzed  the  effects  of  film  mass  and  thickness  in  multilayered  structures 
and  found  strategies  for  maximizing  and  minimizing  phonon  transmission  at  interfaces  and  in  superlattices  [10]. 

Experiments:  effects  of  interfacial  films  on  thermal  conductance 

•  Measured  the  effective  conductance  Hbd  of  Pt/Ni/Si  samples  with  varying  Ni  thickness.  Maximum  enhancement 
of  Hbd  occurred  with  the  thinnest  film  (5  nm),  145%  higher  than  in  the  baseline  Pt/Si  interface:  ^235  vs. 
95  MW  m“^  K“^.  The  dramatic  enhancement  is  partially  due  to  bonding,  but  also  partially  due  to  vibrational 
matching  as  predicted  by  simulation.  Results  to  be  published. 

•  Measured  Pt/Ni/Ge  samples  as  a  counterexample.  In  the  best  case  (5  nm  Ni),  the  effective  Hbd  is  only  50% 
higher  than  that  of  the  baseline  Pt/Ge  interface:  ^160  vs.  105  MWm“^  K“^.  The  smaller  enhancement  is 
consistent  with  an  increase  due  to  Ni/Ge  bonding  in  competition  with  a  decrease  due  to  vibrational  mA-match. 
Results  to  be  published. 

The  grant  also  supported  several  findings  related  tangentially  to  the  proposed  work. 

•  Investigated  thermal  transport  in  systems  involving  carbon  materials.  Determined  scaling  laws  for  conductivity 
of  carbon  nanotube  networks  [11].  Modified  the  DMM  to  predict  Hbd  at  metal-graphite  interfaces  [12],  and 
showed  experimentally  that  surface  treatment  can  cause  Hbd  of  Au-graphite  interfaces  to  vary  by  3  x  [13,  14]. 

•  Simulated  the  vibrational  contribution  to  the  thermal  conductivity  of  Si  near  Tmeit-  Vibrational  energy  transport 
contributes  <  4%  of  the  total  conductivity  of  liquid  Si,  with  electronic  transport  contributing  the  rest  [15,  16]. 

•  Simulated  the  laser  irradiation  of  a  thin  Ag  film  on  Cu  substrate,  leading  to  superfast  melting/recrystallization  of 
the  subsurface  Ag/Cu  interface.  Observed  the  formation  of  a  novel,  lattice-mismatched  interfacial  microstruc¬ 
ture,  effectively  hardening  the  interface  by  laser  irradiation  [17]. 

•  Simulated  the  effect  of  long-range  order  parameter  on  the  thermal  conductivity  of  crystalline  substitutional 
solid  solutions.  Found  that  the  long-range  order  parameter  is  an  effective  method  for  controlling  the  thermal 
conductivity  in  crystalline  substitutional  solid  solutions  (e.g.,  Si^^Gei-^^  alloys)  [18]. 

•  Reviewed  the  research  literature  regarding  methods  for  tuning  heat  transport  at  solid-solid  interfaces,  superlat¬ 
tices,  and  heterostructures  [19]. 

These  findings  have  been  published  in  the  16  archival  journal  papers  and  3  conference  proceeding  papers  cited  here, 

not  including  publications  still  under  preparation.  Each  of  these  publications  acknowledged  funding  from  the  AEOSR. 
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1  Introduction 


The  objective  of  this  research  was  to  investigate  how  the  physical  aspects  of  a  solid-solid  interface  affect  its  thermal 
conductance  /ibd,  which  is  expressed  in  units  of  Researchers  also  often  refer  to  thermal  boundary 

resistance,  which  is  the  inverse  quantity.  As  sketched  in  Figure  1,  /ibd  represents  an  interface’s  capability  to  transmit 
heat  flux  under  a  given  temperature  difference.  Typical  /ibd  at  phonon-mediated  interfaces  (i.e.,  not  metal-metal)  is  in 
the  range  of  10^-10^  MW  m“^  K“^. 


Figure  1:  Sketch  of  a  steady-state  temperature  profile,  T{x),  in  response  to  a  constant  heat  flux,  q,  normal  to  the 
interface  between  two  materials.  The  temperature  drop  at  the  interface,  AT,  defines  its  thermal  boundary  conductance, 
/ibd-  The  inverse  quantity  is  known  as  the  thermal  boundary  resistance  (TBR). 


Understanding  h^D  is  crucial  in  practical  applications.  For  example,  the  GaN-substrate  conductance  is  a  limiting 
factor  in  the  development  of  GaN-based  high-power,  high-frequency  devices  of  interest  to  the  AFOSR  [20-23].  Un¬ 
derstanding  Hbd  has  also  proven  to  be  a  very  interesting  (i.e.,  difficult)  problem  in  fundamental  science,  as  evidenced 
by  active  interest  in  the  topic  from  the  1940s  to  the  current  day  [19].  This  grant  supported  complementary  compu¬ 
tational,  theoretical,  and  experimental  investigations  of  methods  for  controlling  Hbd  at  an  interface  by  modifying  its 
properties.  The  methods,  results,  and  conclusions  of  each  investigation  are  summarized  in  the  following  sections. 


2  Simulations:  effects  of  interfacial  films  on  thermal  conductance 

The  hypothesis  that  motivates  this  grant — that  material  can  be  added  at  an  interface  to  enhance  its  effective  thermal 
conductance — was  first  explored  using  simulations.  Our  primary  tool  was  the  classical  molecular  dynamics  (MD) 
method.  The  first  part  of  the  simulation  work  demonstrated  that  an  interfacial  film  can  enhance  conductance  in  simple 
systems.  The  second  part  laid  the  groundwork  to  extend  those  simulations  to  more  complex  material  systems. 

2.1  Enhancing  conductance  between  simple  fee  solids 

For  details,  see  Refs.  [1]  and  [2]. 

Methods 

This  work  used  molecular  dynamics  (MD)  simulations  to  predict  the  effects  of  inserting  a  thin  film  at  an  interface 
between  Ar  and  “heavy  Ar.”  For  each  configuration,  the  effective  thermal  conductance  Hbd  was  calculated  using  the 
non-equilibrium  molecular  dynamics  (NEMD)  method.  We  prescribed  a  constant  heat  flux  (0.285  GW  m“^)  through 
the  simulation  domain,  which  leads  to  the  development  of  a  steady-state  temperature  profile,  from  which  Hbd  is 
readily  calculated  as  sketched  in  Figure  1 .  The  conductance  was  calculated  in  this  manner  for  a  range  of  masses  and 
thicknesses  of  the  film.  Once  conductance  trends  were  identified,  we  calculated  the  local  density  of  states  in  each 
material  (Ar,  heavy  Ar,  and  the  film),  which  provided  qualitative  explanations  of  the  conductance  trends  based  on 
phonon  transport. 
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Results  and  Conclusions 

A  large  set  of  conductance  calculations  is  presented  in  Figure  2,  each  normalized  relative  to  the  conductance  of  the 
baseline  Ar-heavy  Ar  interface  with  no  film.  From  these  data  we  drew  several  conclusions: 

•  With  interatomic  forces  the  same  throughout,  the  highest  enhancement  occurs  when  the  film  mass  is  roughly 
average.  However,  the  enhancement  is  not  extremely  sensitive  to  mass  near  the  optimum  value,  which  bodes 
well  for  materials  selection.  In  real  systems  with  inhomogeneous  bonding,  we  expect  that  this  trend  with  mass 
corresponds  to  a  more  general  trend  with  vibrational  impedance. 

•  The  enhancement  decreases  monotonically  with  film  thickness — for  maximum  /ibd,  make  the  film  not  thicker 
than  a  few  nanometers. 

•  Mixing  at  the  interface  can  increase  or  decrease  Hbd',  there  exists  an  “optimum”  degree  of  mixing,  which  is 
shown  in  the  right  column  of  Figure  2.  The  additional  enhancement  from  mixing  is  on  the  same  order  as  the 
enhancement  from  the  interfacial  film  itself. 

•  The  enhancement  due  to  the  interfacial  film  is  significant  at  T  =  10%  Tmeit  but  vanishes  at  high  temperatures. 
We  attribute  this  to  shorter  phonon  mean  free  paths  as  temperature  increases,  which  changes  the  transport  in  the 
interfacial  film  from  ballistic  to  diffusive. 


Sharp  interface  Optimum  mixing 


.  .j 

!<'?888888888i 

Interfacial  Film  Mass  (amu) 


Figure  2:  Thermal  conductance  of  various  interfaces  with  interfacial  films,  normalized  by  the  conductance  of  the 
baseline  Ar-heavy  Ar  interface.  The  data  span  four  independent  variables:  film  mass,  film  thickness,  degree  of 
interfacial  mixing,  and  overall  system  temperature. 
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2.2  Transition  to  realistic  interfaces:  metals  and  semiconductors 

For  details  about  the  formulation  of  the  angular-dependent  embedded  atom  method  (A-EAM),  and  for  parameters 
describing  the  binary  Au-Si  and  Au-Ge  systems,  see  Ref.  [3].  For  a  parameterization  of  the  Al-Si  system,  also  see 
Ref.  [24],  which  was  not  supported  by  this  grant.  We  also  present  a  new  preliminary  parameter  set  for  the  Al-Si 
system,  not  yet  published,  which  can  be  used  with  the  Au-Si  parameter  set  to  describe  the  ternary  Au-Al-Si  system. 

Methods 

The  general  form  of  the  A-EAM  was  developed  by  combining  the  existing  Stillinger-Weber  (SW)  potential  for  tetrahe- 
drally  coordinated  materials  [25]  and  the  embedded-atom  method  (EAM)  potential  for  metals  [26,  27].  The  resulting 
potential  functions  are  models  for  interatomic  forces  among  both  covalent  and  metallic  species;  for  example,  in  the 
binary  Au-Si  system.  (A  potential  must  be  available  for  each  pure  component,  either  in  the  SW  form  for  a  covalent 
material  or  in  the  EAM  form  for  a  metal.)  The  new  A-EAM  potential  requires  the  appropriate  selection  of  values  for 
14  parameters  for  each  binary  system. 

For  each  binary  system,  those  parameters  must  be  chosen  to  reproduce  reference  data  of  interest  to  the  simulator. 
In  our  case,  we  are  interested  in  reproducing  interatomic  interactions  in  the  solid  state.  Unfortunately,  little  experi¬ 
mental  information  is  available  for  systems  like  Au-Si  that  are  immiscible  in  the  solid  state,  so  we  performed  sets  of 
calculations  using  density  functional  theory  (DFT)  to  supplement  experimental  data.  The  calculations  in  Ref.  [3]  were 
performed  using  the  Vienna  ab  initio  simulation  package  (VASP),  and  the  calculations  for  the  ternary  Au-Al-Si  system 
(presented  below,  not  yet  published)  were  performed  using  Quantum  Espresso.  The  latter  calculations  were  converged 
within  1  X  10“^  Ryd  with  respect  to  wave  function  cutoff  energy,  energy  density  cutoff,  and  /c-point  sampling. 

The  A-EAM  parameters  were  then  chosen  to  best  reproduce  the  lattice  constants  a,  cohesive  energies  E,  and  bulk 
moduli  B  of  crystalline  structures  as  well  as  formation  energies  of  vacancies  E^vac  and  substitutional  impurities  E^sub- 
In  practice,  three  of  the  parameters  are  held  constant,  as  described  in  Ref.  [3],  so  each  binary  system  has  1 1  free  fitting 
parameters. 

Results  and  Conclusions 

•  A  new  interatomic  potential,  the  A-EAM,  was  developed  for  multicomponent  systems  of  mixed  metallic  and 
covalent  species. 

•  Using  a  series  of  DET  calculations  as  fitting  targets,  parameters  have  been  determined  for  interactions  between 
Au-Si,  Au-Ge,  and  Al-Si  atoms,  as  summarized  in  Table  1.  Those  potential  parameters  give  rise  to  simulated 
material  properties  listed  in  Table  2,  where  they  are  compared  with  reference  values  (either  experimental  or 
DET).  The  properties  of  three  additional  structures,  omitted  here  for  brevity,  were  also  used  for  the  fitting 
process:  AuSi,  AlSi,  and  AU4AI.  In  addition  to  reproducing  behavior  in  the  solid  state,  the  potentials  are  also 
shown  to  accurately  reproduce  enthalpies  of  mixing  in  the  liquid  state. 

•  The  Al-Si  parameters  reported  here  are  compatible  with  the  Au-Si  parameters  and  can  be  used  simultaneously 
to  describe  the  ternary  Au-Al-Si  system  by  applying  the  alloying  rules  described  by  Zhou  et  al  [28].  This  is 
not  true  of  the  Al-Si  parameters  presented,  for  example,  in  Ref.  [24]. 


3  Theory:  understanding  conductance  via  phonon  transport  theory 

To  investigate  the  phonon  transport  underlying  the  conductance  trends  observed  in  Section  2,  we  used  various  theo¬ 
retical  approaches  to  understand  fundamental  phonon  transport  in  interfacial  structures:  (1)  semi-empirical  methods 
such  as  the  DMM,  (2)  the  wavelet  transform  applied  to  MD  simulations,  (3)  Green’s  functions,  and  (4)  the  interfering 
particle  model  (IPM). 

3.1  Modifying  the  diffuse  mismatch  model:  single  interfaces 

Regarding  the  use  of  full  phonon  dispersions  instead  of  the  Debye  approximation,  see  Ref.  [5].  Eor  details  about 
modifying  the  diffuse  mismatch  model  (DMM)  to  include  the  effects  of  optical  phonons,  see  Ref.  [4].  Regarding  the 
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Table  1 :  A-EAM  Parameters  and  Fitted  Values  for  Selected  Binary  Systems 


Parameter 

Au-Si 

Au-Ge 

Al-Si 

Ac  (eV) 

5.189 

4.769 

3.414 

Be  (eV) 

9.8129 

8.779 

5.751 

ac 

-0.555 

-0.555 

-0.517 

Pc 

-0.93 

-0.93 

-0.637 

cre  (A) 

1.81 

1.884 

2.313 

Pc 

4 

4 

4 

Qc 

0 

0 

0 

ymetal 

4.8 

4.2 

1.120 

ycovalent 

16 

16 

16 

^,metal 

ic 

3.56 

3.56 

3.499 

^,covalent 

JC 

-1.1 

-1.1 

-3.130 

ometal 

lie 

2.47 

2.57 

3.229 

Rs  (A) 

2.97 

3.07 

3.073 

Rm  (A) 

4.7837 

4.98 

4.390 

Table  2:  Evaluation  of  Selected  Material  Properties  Predicted  by  the  A-EAM 


Binary  System 

Structure 

Property 

A-EAM 

Reference 

Au-Si 

AuaSi 

a  (A) 

4.193 

4.145^ 

[3] 

4.087^ 

E(cW) 

-4.035 

-3.149* 

-4.108* 

B  (Mbar) 

1.591 

1.203* 

1.888* 

Al-Si 

AlsSi 

a 

3.848 

4.145* 

4.087* 

E 

-4.252 

-3.149* 

-4.108* 

B 

1.172 

1.203* 

1.888* 

Au-Al 

AuAl2 

a 

6.102 

6.070* 

[29] 

E 

-3.416 

-3.149* 

[30] 

B 

0.794 

1.030* 

[29] 

*  Experimental 
^  DFT  calculation  using  EDA 
^  DFT  calculation  using  GGA 
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relaxation  of  the  assumption  of  detailed  balance,  see  Ref.  [6] .  Regarding  the  effects  of  inelastic  phonon  transmission, 
see  Ref.  [7]. 

Methods 

The  DMM  predicts  the  conductance  at  an  interface  based  on  the  carrier  energy  huj,  density  of  states  D,  occupation  /o, 
group  velocity  v,  and  transmission  probability  r  of  each  phonon  frequency  [31]: 

^  X!  ^  dw- 

The  problem  then  becomes  how  to  calculate  r,  since  the  other  quantities  are  based  only  on  bulk  material  properties. 

Results  and  Conclusions 

•  Using  full  phonon  dispersions  rather  than  the  usual  Debye  approximation  can  lead  to  such  and  such  improvement 
in  predictions  of  h^D  in  some  systems. 

•  Contrary  to  common  assumption,  the  transmission  of  optical  phonons  contributes  a  large  fraction  toward  the  total 
conductance  at  many  interfaces  (Figure  3,  left).  The  optical  phonon  contribution  correlates  with  the  number  of 
optical  branches  in  the  abutting  materials;  i.e.,  with  the  complexity  of  their  unit  cells. 

•  The  assumption  of  strict  detailed  balance  when  calculating  the  transmission  probability  r  can  lead  to  poor 
predictions  of  /ibd-  Predictions  can  be  significantly  improved  by  assuming  “quasi-elastic”  scattering  (Figure  3, 
right). 

•  We  developed  an  anharmonic  inelastic  model  (AIM)  which  modifies  the  DMM  to  account  inelastic  4-phonon, 
5-phonon,  etc.  processes  in  a  detailed  manner.  Agreement  with  experimental  /ibd(^)  data  is  significantly  im¬ 
proved. 
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Figure  3:  Findings  regarding  the  phonon  transport  underlying  /ibd,  based  on  the  DMM.  Left:  the  fraction  of  total  Hbd 
contributed  by  the  transmission  of  optical  phonons  at  various  interfaces,  sorted  by  the  total  number  of  optical  branches 
in  the  constituent  materials.  Right:  comparison  of  experimental  /ibd(^)  data  with  DMM  predictions  using  different 
assumptions  for  calculating 


3.2  The  wavelet  transform:  single  interfaces 

For  details,  see  Ref.  [9]. 


7 


Methods 


The  wavelet  transform  is  akin  to  the  Fourier  transform,  except  that  the  basis  function  used  for  the  transform  is  localized 
in  both  position  and  wavenumber.  Thus,  the  wavelet  transform  can  reveal  the  phonon  dynamics  that  are  implicitly 
present  in  a  molecular  dynamics  simulation.  The  wave-packet  method  was  used  to  study  anharmonic  phonon  decay 
and  one-dimensional  NEMD  was  used  to  study  phonon  dynamics  during  transport  at  an  interface. 

Results  and  Conclusions 

The  wavelet  transform  makes  no  assumptions  regarding  the  underlying  physics  of  the  simulation,  and  thus  accommo¬ 
dates  systems  where  anharmonicity  is  important.  Figure  4  shows  the  application  of  the  wavelet  transform  to  analyze 
phonon  dynamics  in  two  scenarios: 

1 .  Left:  the  anharmonic  decay  of  a  large- amplitude  phonon.  The  wavelet  transform  reveals  how  the  energy  distri¬ 
bution  among  wavenumbers  changes  in  time  and  space.  Energy  first  transfers  into  the  2kp  range,  representing 
a  3-phonon  normal  scattering  process.  More  chaotic  interactions  then  occur,  distorting  the  initial  wave-packet’s 
shape  and  introducing  energy  at  greater  wavenumbers.  The  wave-packet  method  is  quite  artificial  and  meant 
to  isolate  a  particular  physical  process.  But  the  wavelet  transform  can  also  accommodate  a  finite-temperature 
system,  which  can  be  thought  of  as  an  ensemble  of  wave-packets. 

2.  Right:  the  distribution  of  energy  among  vibrational  modes  near  an  interface  in  thermal  non-equilibrium.  The 
interatomic  forces  in  the  system  are  purely  harmonic.  The  wavelet  transform  indicates  that  the  cause  of  thermal 
boundary  resistance  is  the  population  of  phonons  in  side  1  that  lie  above  the  cutoff  frequency  of  side  2,  corre¬ 
sponding  to  g  ^  O.dgmax-  Instead  of  transmitting,  these  phonons  accumulate  in  side  1.  The  wavelet  transform 
thus  illustrates  the  underlying  physics  of  the  boundary  resistance,  which  this  grant  has  addressed. 
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Figure  4:  Applications  of  the  wavelet  transform  to  discern  phonon  dynamics  during  MD  simulations.  Left:  the 
anharmonic  decay  of  vibrational  energy  in  a  wave  packet  into  modes  of  different  wavenumbers.  Right:  the  spatial 
variation  of  the  vibrational  energy  spectrum  during  a  NEMD  simulation  at  an  interface. 
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In  conclusion: 


•  The  wavelet  transform  gives  spatially  localized  spectral  information  from  molecular  dynamics  simulations. 

•  It  is  a  flexible  tool  that  can  be  used  to  investigate  individual  scattering  events  (phonon-phonon,  phonon-interface), 
finite-temperature  phonon  non-equilibrium,  and  even  carrier  dynamics  in  non-crystalline  materials. 

3.3  Green’s  functions:  interfaces  with  impurities 

For  details,  see  Ref.  [8]. 

Methods 

The  formulation  for  Green’s  functions,  as  applied  here,  is  based  on  a  system  of  classical  harmonic  Newtonian  equations 
which  describes  the  forces  acting  on  each  atom  in  a  one-dimensional  system.  The  model  system  is  comprised  of  three 
parts:  the  source,  the  sink,  and  a  center  region  (the  “channel”).  In  general,  the  source  and  sink  are  two  different  semi¬ 
infinite  materials,  having  different  masses  and  spring  constants,  and  the  channel  is  some  finite  sequence  of  masses  and 
springs. 

Using  the  Green’s  function  method  and  the  Caroli  formula,  we  solve  the  system  of  Newton’s  equations  to  produce 
analytical  expressions  for  the  energy  transmission  ratio  as  a  function  of  vibrational  frequency.  Those  expressions  were 
used  to  study  the  effects  of  tuning  “impurity”  masses  and  springs  on  phonon  transmission. 

Results  and  Conclusions 

•  At  an  interface  between  materials  1  &  2  with  an  impurity  mass  mint,  the  energy  transmission  of  all  phonon 
frequencies  is  maximized  when  mint  =  (mi  +  m2)/2,  regardless  of  ki  and  k2  (Figure  5,  left). 

•  In  that  case,  the  analytical  expression  for  the  transmission  takes  precisely  the  form  of  the  acoustic  mismatch 
model  (AMM),  which  predicts  transmission  based  on  the  acoustic  impedance  of  each  material.  However,  this 
Green’s  function  approach  leads  to  a  different  expression  for  impedance,  leading  to  a  transmission  expression 
that  is  true  for  all  phonon  frequencies,  unlike  the  AMM. 

•  With  an  impurity  spring  /cint,  the  energy  transmission  of  all  phonon  frequencies  is  maximized  when  ki^t  = 
2  ki  k2/{ki  +  fe)  (Figure  5,  right). 

•  In  that  case,  the  analytical  expression  for  the  transmission  does  not  take  the  form  of  the  AMM.  Rather,  the 
expression  includes  an  additional  resonance  term  that  depends  on  the  properties  of  the  source  and  sink. 


Figure  5:  Phonon  transmission  spectra,  with  curves  corresponding  to  different  impurity  masses,  mint  (left  plot)  or 
different  bonding  strengths,  /cint  (right  plot). 
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3.4  The  interfering  particle  model:  interfaces  with  intermediate  layers 

For  details,  see  Ref.  [10]. 

Methods 

The  interfering  particle  model  (IPM)  [32]  was  used  to  simulate  the  propagation  of  vibrational  energy  through  interfa¬ 
cial  structures  with  N  intermediate  layers.  In  the  original  implementation  of  the  IPM,  Schelling  and  Phillpot  obtained 
the  phonon  transmissivity  at  each  of  the  N  —  1  interfaces  by  using  MD  wave  packet  simulations,  which  were  by  far 
the  dominant  computational  expense.  In  the  current  work,  the  phonon  transmissivities  were  obtained  instead  from 
an  algebraic  expression  that  accurately  reproduces  the  MD  results  but  with  negligible  computational  expense.  This 
permitted  a  large  parameter  sweep  of  layer  masses  and  phonon  frequencies  in  structures  with  =  1  (93  simulations), 
N  =  2  (961  simulations),  and  =  4  (923  521  simulations).  Each  simulation  provides  the  net  energy  transmission 
ratio  for  a  particular  phonon  frequency  through  a  particular  multilayer  structure.  The  results  were  searched  to  identify 
effective  strategies  for  maximizing  and  minimizing  phonon  transport  through  mulitayer  structures. 

Results  and  Conclusions 

Selected  transmissivity  results  are  shown  in  Figure  6  in  structures  with  N  =  1  intermediate  layer  (left)  and  structures 
with  N  =  4  layers  (right).  The  transmissivity  in  the  N  =  4  structure  is  calculated  for  the  low-frequency  case, 
uj  =  4.59  Trads“^.  From  these  data  we  drew  several  conclusions. 


Mi=Ms-  M2  M3=4M, 


Description 

M2 

Ms 

M4 

Ms 

Trans. 

smooth,  geometric 

1.3 

1.7 

2.3 

3.0 

0.976 

smooth,  linear 

1.6 

2.2 

2.8 

3.4 

0.973 

smooth,  opp.  concavity 

2.0 

2.7 

3.3 

3.7 

0.961 

single  interface 

1.0 

1.0 

4.0 

4.0 

0.888 

zigzag,  max.  mismatch 

4.0 

1.0 

4.0 

1.0 

0.654 

zigzag,  min.  trans. 

3.9 

1.0 

4.0 

1.0 

0.623 

Figure  6:  Selected  phonon  transmissivity  data  predicted  using  the  IPM  in  structures  with  A^  =  1  intermediate  layer 
(left  plot)  and  in  structures  with  A^  =  4  layers  (right  table).  In  both  cases,  the  left  lead  has  normalized  mass  Mi  =  1, 
and  the  right  lead  has  M7V+2  =  4. 


•  In  the  structure  with  N  =  1  layer,  choosing  the  layer  mass  as  the  geometric  mean  maximizes  the  energy 
transmission  of  low-frequency  phonons.  For  higher  phonon  frequencies,  the  optimum  choice  of  mass  shifts 
upward  toward  the  arithmetic  mean.  Therefore,  the  choice  of  layer  mass  that  maximizes  the  effective  /ibd  is 
likely  in  between. 

•  In  the  structure  with  A^  =  4  layers,  transmission  is  maximized  when  the  layer  masses  change  monotonically 
(labeled  “smooth”).  The  optimum  transmissivity  occurs  when  the  sequence  of  masses  is  geometric,  providing 
a  significant  enhancement  compared  to  the  single-interface  transmission.  However,  this  enhancement  is  not 
particularly  sensitive  to  the  exact  sequence  as  long  as  the  transition  is  “smooth.” 
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•  The  transmission  is  minimized  when  the  layer  masses  alternately  increase  and  decrease  (“zigzag”),  as  in  a 
superlattice.  However,  the  minimum  transmissivity  does  not  necessarily  occur  in  the  sequence  with  the  greatest 
impedance  mismatch  (“max.  mismatch”),  since  some  constructive  interference  can  occur.  Instead,  the  minimum 
transmissivity  occurs  in  an  aperiodic  structure  which  sacrifices  some  degree  of  impedance  mismatch  to  disrupt 
constructive  interference. 


4  Experiments:  effects  of  interfacial  films  on  thermal  conductance 

Finally,  the  findings  from  simulations  and  theoretical  analysis  were  used  to  design  a  series  of  experimental  measure¬ 
ments  of  /ibd  at  interfaces  with  varying  thicknesses  of  interfacial  films. 


Methods 

The  effective  thermal  boundary  conductance  of  each  metal-substrate  system  was  measured  using  time  domain  ther- 
morefiectance  (TDTR),  an  optical  pump-probe  technique  that  is  well  established  for  measuring  thermophysical  prop¬ 
erties  of  systems  with  nanoscale  features  [33].  Two  series  of  samples  were  fabricated  at  the  Army  Research  Laboratory 
in  collaboration  with  the  US  Naval  Academy:  (1)  a  series  of  Pt-Ni-Si  samples,  in  which  the  thickness  of  the  Ni  layer 
was  varied  from  5  to  100  nm,  and  (2)  a  corresponding  series  of  Pt-Ni-Ge  samples.  The  choice  of  materials  was  based 
in  part  on  the  ratios  of  the  Debye  temperatures  Od  of  the  intermediate  layer  and  the  substrate  [2,  34].  The  thermal 
decay  of  each  sample,  measured  by  TDTR,  was  compared  to  a  multilayer  solution  of  the  heat  diffusion  equation  to 
extract  the  conductance  /ibd  between  Ni  and  Si(Ge)  [35].  Using  a  series  thermal  resistor  analysis,  the  total  Pt-to- 
substrate  conductance  is  calculated.  That  total  conductance  is  compared  with  the  baseline  measurements  of  /ibd  at 
Pt-Si  and  Pt-Ge  interfaces  without  an  intermediate  layer. 


Results  and  Conclusions 

Figure  7  shows  the  resulting  /ibd  measurements  as  a  function  of  intermediate  layer  thickness  for  each  of  the  two 
material  systems.  From  these  results  we  drew  several  conclusions: 

•  The  Pt-Ni-Si  system  shows  drastically  enhanced  conductance  compared  to  the  Pt-Si  interface.  The  enhance¬ 
ment  is  best  for  the  5  nm  Ni  layer  and  decreases  with  increasing  thickness,  as  predicted  by  the  MD  simulations 
summarized  in  Section  2. 

•  Some  of  the  enhanced  /ibd  can  be  attributed  to  increased  bonding  between  Ni-Si  (vs.  Pt-Si).  The  remainder  of 
the  enhancement  is  due  to  the  increased  vibrational  match:  the  Ni  layer  has  a  Debye  temperature  in  between  the 
Pt  and  Si  leads. 

•  In  the  Pt-Ni-Ge  system,  /ibd  is  enhanced  slightly  when  the  Ni  layer  is  very  thin.  However,  as  thickness  in¬ 
creases,  quickly  falls  below  the  baseline  Pt-Ge  conductance.  We  attribute  this  behavior  to  enhanced  Ni-Ge 
bonding  in  competition  with  a  vibrational  mismatch:  the  intermediate  Ni  layer  has  a  higher  Debye  temperature 
than  both  of  the  Pt  and  Ge  leads. 

•  The  large  enhancement  observed  in  the  Pt-Ni-Si  system  and  the  small  enhancement  in  the  Pt-Ni-Ge  system 
provide  strong  evidence  that  /ibd  can  be  enhanced  at  mismatched  interfaces  by  inserting  additional  material  of 
intermediate  properties. 

The  experimental  work  summarized  here  is  under  preparation  to  be  published. 
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Figure  7:  Effective  conductance  of  Pt-Ni-semiconductor  samples  with  varying  thickness  of  the  Ni  layer.  The  Ni  layer 
is  expected  to  “bridge”  the  vibrational  mismatch  betwen  Pt  and  Si  (left),  but  not  between  Pt  and  Ge  (right). 
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